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Abstract 

We derive the explicit analytic expression for the Stokes' drift in one 
dimension in the presence of a dichotomic Markov forcing. For small 
amplitudes of the forcing, the drift is enhanced, but the enhancement is 
reduced with increasing frequency of the forcing. On the other hand, a 
reduction of the drift or even a flux reversal can be induced at larger 
amplitudes, while the flux is now found to be an increasing function of 
the perturbation frequency. 

PACS numbers: 02.50.-r, 05.40.+j, 05.60,+w 

1 Introduction 

A longitudinal wave travelling through a fluid imparts a net drift motion to 
the suspended particles - an effect known as Stokes' drift. The classical Stokes' 
drift 0| refers only to the deterministic behavior, i.e. it does not account for the 
stochastic fluctuations or perturbations in the system. It has been extensively 
studied in various practical contexts like, for instance, the motion of tracers in 
meteorology and oceanography Q and that of the (doping) impurities in crystal 
growth [||. A simple intuitive explanation of classical Stokes' drift is that the 
suspended particles spend a longer time in the regions of the wave-train where 
the force due to the wave acts in the direction of wave's propagation than 
the time spent in those regions where the force acts in the opposite direction; 
therefore particles are driven in the sense of wave's propagation. When several 
linearly superposed waves are present, the resulting drift velocity is simply the 
sum of the contributions from each wave . 

Recently, there has been some interest in stochastic Stokes' drift (|^, ||, |j, see 
also [§). It was found that the diffusive motion of the suspended particles can 
significantly alter the amplitude of the drift velocity. But it is the merit of H 
to open stochastic Stokes' drift to further investigation by putting it in its nat- 
ural context, namely ratchet-like Brownian motors in which the thermal motion 
of small (microscopic) particles is rectified by an asymmetric time-dependent 
potential This scenario is referred to as the flashing ratchet. Paper Q is 
limited to the case when the motion of the overdamped suspended particles is 
diffusion-dominated, i.e. when the deterministic forcing due to the wave is very 
weak. An approximate perturbative technique is developed for calculating the 
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drift velocity in arbitrary dimensions, which shows that the drift depends on 
particles' diffusivity in magnitude and, in more than ID, also in direction (a very 
appealing feature, useful in particle selection devices - see, for example, in fl 
Doering et all, 1994, and in || M. Bier, 1996). Papers Q and, in a very different 
context Q|, give an exact integral expression for Stokes' drift of overdamped dif- 
fusive particles in ID, for arbitrary wave forms. Using a Fokker-Planck equation 
technique, it finds that white noise reduces Stokes' drift by comparison to its 
classical (deterministic) value, the flux reversal being absent. 

The purpose of this short communication is to show that the characteristics 
of the ID Stokes' drift may become more complex (including flux reversal or 
enhancement of the classical drift) when the particles are subjected to stochastic 
forcing that is non- white (i.e. it has a finite correlation time). We will focus on 
the case of dichotomic Markov forcing, in which the particle is subjected to a 
randomly alternating force, in addition to the force exerted by the wave. For this 
choice of coloured noise, explicit analytic results can be obtained, while this is 
not the case for other types of coloured noise, including the Ornstein-Uhlenbeck 
process. Even though the time-average of the dichotomic force is taken to be 
zero, its effect in conjunction with the asymmetry of the wave's perturbation 
is highly nontrivial. As such the model realizes in the Stokes' drift context 
another paradigm of Brownian motors, namely that of the so-called rocking 
ratchet. In this scenario, the systematic motion acquired under influence of 
an alternating zero-average force essentially revolves on the asymmetry of the 
nonlinear response in the presence of a steady but asymmetric potential ||. 
The dichotomic Markov forcing presents another important advantage in that 
the exact integral expression for Stokes' drift can be derived. Furthermore, in 
some cases of interest, the resulting integrals can be evaluated explicitly. Our 
analytic results are also complemented with computer simulations, which are 
found to be in full agreement with the theory. 



2 The model 

The starting point is the following stochastic differential equation: 

x t = f(x t -vt)+at), (i) 

where £(<) is a symmetric dichotomic noise (a two-step Markov process) |Tcj . It 
can take only two values, ±A , with equal probability, jumping between them 
with a probability k per unit time. It has zero mean and its autocorrelation 
function (£(t) £(t')) — A 2 exp(— 2 k\t — t'\) shows a finite correlation time r c = 
l/2k. 

With a suitable time scaling, eq.([lj) models the ID overdamped motion of a 
Brownian particle subjected to a dichotomic forcing. X t is the particle position 
at time t, while / is the periodic forcing due to the wave travelling at the speed 
v with wavelength A: 

f(y + X) = f(y),Vy, (2) 

with f(y) < v , Vy as the case of physical relevance. 

The quantity of interest is the drift velocity, which represents the mean 
velocity of the particles in the long time limit: 

= lim ~ (X t ) , (3) 

t— >oo t 
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where the mean is to be taken over the realisations of the dichotomic noise. By 
defining a new variable Y t = X t — vt one can rewrite eq.(|l|) as: 

Y t = F(Y t )+at) , (4) 
with F(y) = f(y) — v, and therefore: 

i? — v = lim - (Y t ) 

t — >00 t 

/+oo 
dyF(y)[p + (y,t)+p-(y,t)}. (5) 
-oo 

Here p±(y, t) represent the probability densities for the particle position y and 
the value of the dichotomic noise equal to ±A, respectively. Their time evolution 
is described by |To[ : 

^ = ~ [(F-A)p_] -*(,_-*.). (6) 

The first term in the r.h.s. of each of these equations corresponds to the deter- 
ministic flow in phase space (between two jumps of the dichotomic noise), while 
the other term corresponds to the jumps of the dichotomic noise between ±A. 
Because of the periodicity ^ of f(y) (F(y)), in order to evaluate (||) we only 
need to study the asymptotic behavior of the reduced probability density: 

oo 

P{y,t)= [p+{y + n\,t)+p_{y + n\,t)}. (7) 

n — — oo 

From its definition, it is obvious that the reduced probability density is subjected 
to periodic boundary conditions: 

P(y + \t) = P(y,t) , Vy,t, (8) 

while the normalization condition reads: 



f P(y,t)dy = l. 
Jo 



(9) 



Considering also the quantity 



p(y,t)= J2 [p+(y + n\,t)-p^(y + n\,t)\ , (10) 



n— — oo 



it follows from (||) that P(y, t) and p(y, t) obey a set of coupled evolution equa- 
tions: 
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The stationary solution P st {y) can be easily found pif : 
i rv+>> TP'( V '\ + 9fr 

P " (2/) = z i, J 12 ^)-^ cxp[ -^ (?/,) + ^ (12) 

which is valid as long as F 2 (y) — A 2 ^ 0, Vy. Here <^(y) is defined by: 

and i? is the normalization factor according to eq.@. Therefore the drift ve- 
locity can be finally expressed as: 

A 

.St 



f(y)P st (y)dy. (14) 

These results are valid for an arbitrary wave form, provided that the obvious 
necessary differentiability and integrability conditions are fulfilled. Note that 
the drift depends nonlinearly on F, so that the contributions of different waves 
are not additive. 



3 The piecewise linear wave 

An explicit analytical expression for P sf (y) can be obtained in the case of a 
piecewise linear wave forcing: 



f(v) = 



-(l-b)v, y € [0,A/2-2e[mod A 

-[(l-b)v + bv(y- A/2 + 2e)/e], y G [A/2 - 2s, A/2[ mod A 
-(l + b)v, y e [A/2, A - 2s[ mod A 

-[(1 + 6)t> - - A + 2s)/s], y € [A — 2s, A[ mod A, 

(15) 

with < b < 1. A square wave corresponds to s — ► 0, while s = A/4 represents 
the triangular-like wave. 

The final expressions for P st (y) are rather lengthy and will not be reproduced 
here. But we performed numerical simulations, which show very good agreement 
with theoretical results, cf. fig. ([!]). Note that these results are valid for < 
A < (1 - b)v or for A > (1 + b)v. 

The analytic results for the square-like wave, corresponding to the limit 
£ — > 0, are of more interest. They read as follows: 

P st {y) = Z- 1 { - - + - — [e- u -e- u 1 



(l-b)v {l-b 2 )v L 1 (1 - b) 2 v 2 - A 2 

for y G [0, A/2[ , and 



1 - e- u ~ v e^v/A-i) r _ v _ v _ v , 2A 2 b 



(16) 



P St (y) = \ -TT^-TT— ~ Tj— [ e ' -« 



(l + &)« (l-6 2 )v J (l + 6) 2 v 2 - A 2 

(17) 
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for y £ [A/2, A[ , where we introduced the dimensionless variables: 

a(l-b) 



U 
V 



2[(l-6) 2 /3 2 - 1] ' 

a(l + b) 
2[(l + 6) 2 /3 2 -l] ' 

2/ " A P = - A (18) 



A 2 ' r ~ A 
The following expression for the drift velocity is obtained: 

^ i- e ^- II 4^(i + e c/ + v - e c/ - e 1/ ) : 



(19) 



where i? c ; = 6 2 u is the classical value of the Stokes' drift velocity. We mention 
some limiting cases of interest. First, the limit A — > or k — > oo (i.e. (3 — > oo 
and/or a — > oo) leads to the classical value of the Stokes' velocity i?/# c z = 1. 
Second, the drift velocity derived in || is recovered in the white noise limit 
A — » oo, fc — » oo, with v4 2 /2fc = Z? finite (/3 — > 0). Finally, the quenched-noise 
limit fe — ► (a — * 0) (i.e. half of the particles, chosen at random, are subjected 
to a constant external forcing +A, while the other half to an external forcing 
—A) results in the following mean velocity: 



Vol /3 2 - 1 



(20) 



This expression reveals the possibility of a flux reversal: the drift velocity is 
negative for /3 < 1, i.e. for sufficiently large amplitude. 

Turning to the general case, we have represented in Fig. || the results of 
eq. (|l^) together with simulations for the square wave for two representative 
sets of parameter values. For small values of the forcing amplitude A one finds 
that the drift is enhanced (i?/# c z > 1), but the enhancement is reduced upon 
increasing the frequency of the perturbation. In fact one can prove that en- 
hancement occurs whenever A < (1 — b)v, by noting that U, V > and making 
use of the inequality e x > 1 + x. On the other hand, a reduction of the drift or 
even a flux reversal (z9 < 0) can be induced at larger amplitudes (for small k, 
according to eq. (|2^)), while the flux is now found to be an increasing function 
of the perturbation frequency. 



4 Concluding remarks 

In conclusion, the rocking ratchet version for Stokes' drift reveals two new fea- 
tures, which are to be attributed to the colour (finite correlation time) of the 
applied force: the drift is enhanced by small amplitude perturbations, while a 
flux reversal occurs at sufficiently large amplitude and large correlation time. 
While our calculations are limited to the case of dichotomic Markov noise, we 
expect that similar results hold for other coloured noises. 

In H it was reported that, in two dimensions, both the direction and ampli- 
tude of the Stokes' drift can change with the intensity of the Brownian motion 
when several waves are present. A similar phenomenon will appear here when 
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one considers, for example, two waves propagating in orthogonal directions X\ 
and X2- A dichotomic forcing £ will induce different drifts along these direc- 
tions depending on its projections £1 = ±Ai and £ 2 = ±A 2 , as well as on its 
transition rate k. In particular, by varying k one is able to induce a significant 
change in the resulting drift direction, as can be seen in fig. [| for the simple 
example of identical square waves. 
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Figure 1: Stationary probability density for the piecewise linear wave, eq. (|l5|), 
with v = 1, b = 0.5, A = 10 and e = 1 (the arrows indicate the boundary points 
of the piecewise linear regions). Dashed lines represent the theoretical results, 
while solid lines correspond to simulations of 10 6 particles evolving for 200 time 
units. The insets show the two realisations of y. Fig. 0(a): A = 0.2, k = 1. 
Fig. 0(b): A = 2, k = 0.1. 
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Figure 2: Drift velocity according to eq. (j3|) for a square wave with the same 
values of the parameters as in Fig. [|, except here e — 0. Symbols represent 
simulations with 10 3 particles measured after 5 x 10 4 time units. The lines 
represent the theoretical value, eq. (p^). Fig.||(a): drift enhancement. Fig.|2|(b): 
flux reversal. 
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Figure 3: The resulting drift velocity (in units d c i] thick arrows) for two iden- 
tical square waves (with the same parameters as in Fig. propagating respec- 
tively along the + x\ and + X2 axes, in the presence of a stochastic forcing with 
A\ = 0.45, A2 — 1.6 (applied along the direction of the thin arrows), for dif- 
ferent transition rates k. The dashed arrow corresponds to the classical limit 
k — > 00. 
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